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I.  INTRODUCTION 


One  basing  concept  for  the  MX  missile  is  to  have  the  missile  movable 
within  long  underground  tunnels.  In  order  to  estimate  the  survivability  of  this 
system,  it  is  necessary  to  estimate  the  rate  at  which  blast-induced  air  shock 
waves  are  attenuated  within  the  tunnel.  This  requires  a knowledge  of  the  tur- 
bulent wall  boundary  layer  induced  by  the  shock  wave. 

A theory  for  the  turbulent  boundary -layer  development  behind  a shock 
wave  moving  with  uniform  velocity  has  been  presented  and  simplifications 
introduced.  2 Effects  of  wall  blowing  (ablation)  were  not  considered.  The  pur- 
pose of  the  present  study  is  to  summarize  the  previously  reported  theory 
(Appendix  A),  to  present  numerical  results  for  shock  propagation  in  air 
(Table  1),  to  provide  approximate  analytic  expressions  (Appendix  B),  and  to 
provide  estimates  for  effects  of  blowing  (Appendix  D).  The  notation  of  the 
previous  works 2 is  used  and  is  summarized  in  the  Symbols  section.  The 
present  theory  is  valid  for  cases  where  the  increase  in  wall-surface  tempera- 
ture is  small  relative  to  the  free -stream  recovery  temperature.  This  approxi- 
mation is  examined  in  Appendix  C. 


*H.  Mirels,  Boundary  Layer  Behind  Shock  or  Thin  Expansion  Wave  Moving 
into  Stationary  Fluid,  TN  3712,  National  Advisory  Committee  for  Aeronautics, 
(May  19  56^. 

2H.  Mirels,  "Shock  Tube  Test  Time  Limitation  Due  To  Turbulent  Wall  Bound- 
ary Layer,  11  AIAA  J.  2^(1)*  84  (1964). 
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II.  THEORY 


In  Refs.  1 and  2,  the  turbulent  boundary  layer  behind  a shock  moving 
with  uniform  speed  is  considered  in  shock -fixed  coordinates.  The  flow  is 
steady  in  this  coordinate  system  and  is  illustrated  in  Fig.  i.  The  wall-surface 
temperature  is  assumed  to  remain  at  its  initial  value.  Symbols  are  defined 
in  the  Symbols  section. 

The  theoretical  development  of  Refs.  1 and  2 is  outlined  in  Appendix  A. 
Numerical  results  for  the  turbulent  boundary  layer  behind  a shock  propagating 
into  ideal  air  and  real  air^’^  are  given  in  Tables  la  and  lb,  respectively. 

Wall  blowing  effects  are  ignored.  An  initial  air  temperature  T^  = 522°R  is 
assumed.  The  results  are  only  weakly  dependent  on  T^.  Measurements  of 
the  turbulent  heat  transfer  behind  a shock  propagating  in  air  have  been  pre- 
sented for  shock  Mach  numbers  in  the  range  3 s M < 8.  These  results  are 
reproduced  in  Fig.  2.  Comparison  with  the  heat  transfer  predictions  in 
Table  lb  shows  good  agreement  (within  a few  percent  of  the  mean  experimental 
value  at  each  shock  Mach  number).  This  agreement  tends  to  confirm  the 
theory  of  Refs.  1 and  2 for  Mach  numbers  in  the  range  Mg  S 0(10)  and  non- 
blowing  walls. 

Turbulent  boundary -layer  properties  of  major  interest  are  local  bound- 
ary thickness,  6;  wall  shear,  r = (|i.6u/6y)  ; heat  transfer  to  the  wall, 

-q  = (k6T/6y)  ; and  the  vertical  velocity  at  the  edge  of  the  boundary  layer,  v . 
*w  w c 


U.  H.  Lewis  and  E.  G.  Burgess,  in.  Charts  of  Normal  Shock  Wave  Properties 
in  Imperfect  Air,  AEDC-TDR -64-43,  Arnold  Engineering  Development  Center, 
Air  Force  Systems  Command  (March  1964). 

*C.  H.  Lewis  and  E.  G.  Burgess,  III,  Charts  of  Normal  Shock  Wave  Properties 
in  Imperfect  Air  (Supplement:  Ma  = 1 to  10),  AEDC-TR-65-196,  Arnold 
Engineering  Development  Center,  Air  Force  Systems  Command  (September 
1965). 

5R.  A.  Hartunian,  A.  L.  Russo,  and  P.  V.  Marrone,  "Boundary  Layer 
Transition  and  Heat  Transfer  in  Shock  Tubes,  " J.  Aerospace  Sci.  27, 

587  (I960). 
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Table  la.  Turbulent  Boundary  Layer  Behind  Shock, 
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are  based  on  Appendix  A. 
0.72,  hw  = ht,  T j = 522°R; 
in  atmospheres,  x is  in  feet. 
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TURBULENT  HEAT  TRANSFER  IN  SHOCK  TUBE 

Fig.  2.  Comparison  of  Present  Theory  With  Experimental 
Data  of  Ref.  5 
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O EXPERIMENT  - GLASS  WALL 
□ EXPERIMENT  - STEEL  WALL 


1 


Approximate  analytic  expressions  for  these  quantities  are  given  in 
Appendix  B.  These  are,  for  W > 1. 02  and  = 522°R 


(p.x)i/56/x=  0.  0370  j— ^ 

' ( [w  + 


W - 1 


(7/3)] 


1/5  _ -2T„(Plx)‘ 


p (u  - U )' 

re\  w el 


= 2. 22  X 


10-3|-W  111/3)1 
(W  - l)2  _ 


(-Ota  x)’ 


(Plx)1/5cf 


/_  u/5_  _ 1^1 

Plx  CH  = p~ (u  - u )(Ti  ~h  ")  = 2/3 

' Ke  w e r w 2Pr 


ta*)1 


(-%> 


u - u 
w e 


= 4.65  X 10' 


(W  - 1)J 
[W  + (7/3)]' 


where  x is  in  feet,  p^  is  in  atmospheres,  and  W = u^/u^  = Pe/Pj  is  the 


velocity  (density)  ratio  across  the  shock.  Other  symbols  are  defined  in  the 
Symbols  section.  Equations  (la)  through  (Id)  agree  with  the  numerical  results 
in  Table  1 to  within  about  10%.  The  quantities  C^,  C^,  and  vg  are  normalized 


by  the  use  of  u - u , which  is  the  free-stream  velocity  in  laboratory  (wall 

W 0 


stationary)  coordinates.  It  is  seen  that  vg  is  negative,  i.e.,  the  vertical 
velocity  at  the  edge  of  the  boundary  layer  is  directed  toward  the  wall.  This 
is  the  result  of  the  aspirating  effect  of  the  wall  in  shock-fixed  coordinates 
(Fig.  1). 


The  theory  of  Refs.  1 and  2 assumes  that  the  departure  of  the  wall- 
surface  temperature,  T , from  its  initial  value,  T^,  is  negligible.  In 
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Appendix  C,  the  variation  of  wall-surface  temperature  with  distance  behind 
the  shock  is  shown  to  equal  (for  strong  shocks) 

0.6  0.3  fe  ‘ ‘)=em.!/V,/5[‘  +0(w)]  <2*> 

pt  X \ 1 / 

where  E is  a function  of  wall  material;  typical  values  (from  Carslaw  and 
Yaeger6)  are  given  in  Table  2.  The  present  theory  is  valid  when  the  change 
in  wall-surface  temperature  is  small  relative  to  the  recovery  temperature 
T . For  strong  shocks,  the  latter  ratio  is  found  from 


w 


- T. 


“OTT 

Pt  X 


E 

0.38 


m1/2w4^5 

s 


[1+0(ir)] 


(2b) 


The  above  equations  neglect  ablation,  which  tends  to  limit  the  wall-surface 
temperature  increase  and  therefore  to  extend  the  region  of  validity  of  the 
present  results.  However,  the  effects  of  blowing  on  the  boundary -layer  prop 
erties  need  to  be  considered.  These  effects  are  treated  in  Appendix  D. 

The  effect  of  blowing  on  boundary -layer  properties  at  a fixed  x location 
can  be  expressed  by 


CH  Cf  f n(  1 4 B') 

<CH>0 


(3a) 


(1  + B^njl 

B7 


+ B') 


(3b) 


H.  S.  Carslaw  and  J.  C.  Yaeger,  Conduction  of  Heat  in  Solids,  Oxford 
University  Press  (1959),  pp.  76,  496. 
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Table  2.  Approximate  Thermal  Properties  of  the  Wall3, 


Type 

Material 

(pCk)‘/2. 

cal 

E ,b 

/i_\0-8/n0*3 

\atm/  \ ft  / 

2 . TTz 

cm  K (sec) 

Metal 

Aluminum 

5. 17  X 10'1 

3.92  X 10-4 

Mild  steel  (0.  1%  c) 

3.  19.  X 10_1 

6.34  X 10'4 

Cast  iron 

3.48  X 10_1 

5.82  X 10'4 

Nonmetal 

Glass  (crown) 

3.67  X 10‘2 

5. 52  X 10'3 

Concrete  (1:2:4) 

3.41  X 10"2 

5.93  X 10'3 

Gas 

Air 

1.34  X 10'4 

1.51  X 10° 

From  Ref.  6 

bE  = 2.024  X 10-4/(pCk)y2 


where  subscript  zero  denotes  the  local  nonblowing  value  [given  by  Eqs.  (la) 
through  (Id)]  and 


B' 


(pv), 


w 


p ju  - u I CT. 
Ke  I w el  H 


hr 

AH 


(3c) 


Here  AH  is  the  difference  between  the  enthalpy  of  the  unheated  wall  material 
and  the  wall  material  at  the  gas/wall  interface  (i.e.,  the  effective  heat  of 
ablation  of  the  wall  material).  The  effect  of  blowing  on  vg  is  given  in  Eq.  (D-10). 

The  present  results  are  applicable  for  a shock  propagating  in  a tube  of 
diameter  d provided  that  26 /d  « 1 (i.e.,  the  boundary  layer  is  thin  relative 
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to  the  tube  radius).  The  present  results  can  be  taken  as  approximately  correct 
for  values  of  26/d  up  to  about  1.  Boundary -layer  closure  (i.e.,  26/d  = 1) 
occurs  at 


= (Pld)1/4/[2(Plx)1/5(6/x)]5/4 

(4a) 

= 65  ( P±d)1/4(1  ± 0.05)  Mg  > 4 

(4b) 

Equation  (4b)  applies  for  nonblowing  cases.  Blowing  causes  boundary -layer 
closure  to  occur  at  smaller  values  of  x.  The  closure  stations  for  the  blowing 
and  nonblowing  cases  are  related  by  [e.g.,  Eq.  (D-ll)] 


Cc  _r  s'  y 

717  'L(i  + B')f  n(  1 + B')  J 


The  flow  downstream  of  xc  is  fully  developed.  The  shear  and  heat  transfer 
in  this  region  are  on  the  order  of  the  shear  and  heat  transfer  at  xc* 

The  present  results  can  be  used  to  estimate  shock  attenuation  for  cases 
where  26 /d  « 1.  The  effect  of  the  boundary  layer  on  the  inviscid  flow  ex- 
ternal to  the  boundary  layer  is  found  from  a distribution  of  volumetric  mass 

7 

sources  defined  by 


• 4 

M = -v  p v 
d re  t 


mass  addition 
volume  -time 


where  vg  is  evaluated  by  use  of  Eqs.  (id)  or  (D-8).  A combination  of  volu- 
metric momentum  and  energy  sources  is  needed  when  the  flow  is  fully 
developed  and  the  wall  is  nonblowing.  Mass  sources  must  be  included  in  the 


H.  Mirels  and  J.  F.  Mullen,  "Small  Perturbation  Theory  For  Shock  Tube 
Attenuation  and  Non  Uniformity,  " Phys.  Fluids  7 (8),  1208  (1964). 
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HI.  CONCLUDING  REMARKS 


The  present  results  are  intended  to  facilitate  estimates  of  the  turbulent 
boundary  layer  behind  shocks  moving  in  air.  The  results  are  expected  to  be 
fairly  accurate  for  nonblowing  walls  and  for  Mfl  S 0(10)  (e.g.,  Fig.  2).  Further 
study  is  required  in  order  to  validate  the  results  for  Mfl  >10  and  a nonblowing 
wall.  Effects  of  blowing  may  be  considered  uncertain  for  all  Mfl.  Further 
comparisons  with  experiment  and  comparisons  with  numerical  results  from 
more  accurate  turbulent  boundary-layer  codes  are  desirable  for  >10  and 
for  wall  blowing. 


APPENDIX  A 


BOUNDARY- LAYER  CHARACTERISTICS 


Characte  e mi -infinite,  zero-pressure -gradient, 

steady-state  tui"  J«try  layer  are  summarized  herein  for  the  case 

of  a moving  waL.  0.  lions  for  the  boundary  layer  behind  a shock,  moving 

with  uniform  speed,  are  then  noted.  These  results  are  based  on  the  theoreti- 
cal development  in  Refs.  1 and  2. 

1.  GENERAL  MOVING  WALL 


Consider  a zero -pressure -gradient  turbulent  boundary  layer  with  a 
moving  wall  in  a coordinate  system  in  which  mean  properties  are  steady 
state  (e.g.,  Fig.  1).  The  wall  velocity  is  denoteduw,  and  the  free-stream 
velocity  is  denoted  u . Following  the  development  in  Ref.  1,  the  boundary- 
layer  thickness,  6 , and  momentum  thickness,  0,  are  related  by 


9/6  =7he  r1  C6[W  - (2W  - 1)  c + (W  - 1)  £Z1  d C 
w •'0  1 + be  - c£ 


(A- la) 


6*/6 1 

TTW  ~ l.w 


1 . 7 he  f 1 C6  rw  - (W  - i)  C]dC 
^ *'0  1 + be  - c£ 


(A- lb) 


where 


b =ir  - 1 

w 


■ft-) 


(A  -2a) 


(A-2b) 


EDINO  PASS  BUNK-NDT  FILMED 


Corresponding  values  of  6*  and  9 are  found  from  Eqs.  (A- la)  and  (A- lb). 
Noting  that  d6/dx  = (4/5)  (6/x),  we  find  the  shear  and  heat  transfer  from 


w 

T 


P u 
re  e 


_ Pr2/3(W  - 1)  _ d©  _ 4 / 9 \ £ 

p u (Ti  - h qw  dx  5 \ 6 / x 


(A-5) 


e e r w 


where  t = (p  8u/9y)w  and -q^  = (kaT/0y)w«  The  vertical  velocity  at  the 
edge  of  the  boundary  layer,  ve>  is  obtained  from 


e _ d8  _ 4 6. 


— = 7ET"  3 

e 


(A- 6) 


Local  shear  and  heat  transfer  coefficients  may  be  expressed  in  the  form 
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The  latter  are  positive  quantities  and  are  referenced  to  "wall  fixed"  quantities. 
(Note  that  | uw  “ ue  I *8  t^le  velocity  of  the  free  stream  relative  to  the  wall. ) 

It  follows  from  Eqs.  (A- 5)  and  (A-7)  that  C^/2  = Pr^3  C^.  With  the  intro- 
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1.  6 /e  /6  \ 6 
i w - 1 1 iTTwix 

0.8  hlb\b 

,3 |w - n r7 x 


(A-7a) 
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duction  of  a Reynolds  number,  Rg  = | u^  - ue| 


Re 1/5  C 


H 


0.0459 

-77* 


0/6 

TT-w 


x/v  , it  follows  that 
e 

P /P  V 

m e/ 


W - 1| 


(A-8) 


In  the  limit  W = 0,  the  above  equations  define  the  boundary  layer  over  a semi- 
infinite flat  plate.  The  case  W > 1 corresponds  to  flow  behind  a shock  wave 
moving  with  uniform  velocity  and  is  considered  in  the  next  section. 

In  Ref.  1,  Eqs.  ( A - 1 ) through  (A-8)  were  derived  on  the  basis  of  an 
ideal  gas,  namely,  h~  T~  p”  * at  constant  pressure.  The  present  form, 
following  Ref.  2,  attempts  to  reduce  the  sensitivity  to  this  assumption  to 
some  extent. 

2.  BOUNDARY  LAYER  BEHIND  SHOCK 

We  consider  a shock  moving  with  uniform  velocity,  u = u , into  air. 

s w 

Undisturbed  air  conditions  are  denoted  by  subscript  1 (Fig.  1).  We  assume 
Yt  = 7/5,  Pr  = 0.  72,  r(o)  = 0.90,  Tw  = T = 522°R,*  and  hw  = hj.  The  latter 
assumption  is  discussed  in  Appendix  C.  Corresponding  values  of  Mg,  W = 
u /u  = p /p  , and  T /T . are  found  from  ideal  or  real  gas  (Refs.  3 and  4) 

W 0 0 W 0 1 


The  solution  is  relatively  insensitive  to  the  assumption  T,  = 522 °R.  For 
*■  ~ Ta),  it  can  be  shown  6 ~ T 0- 2^  + 0.1. 


1 
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shock  relations.  Noting  0.  7 < (h  /h  ) < 1.0,  we  further  assume  that  T /T.  = 

me  ml 

(h  /h  ) (T  /T.)  and  p / p = h /h  , where  h /h  is  found  from  Eq.  (A- 4) 

' m e e 1 Km  Ke  e m me 

and 


1 + 0.  2 M 


(A-9a) 


1 + 


0.  18  M2 


For  air  at  Tj  = 522®  R 


6.93  X 106 


(atm- ft) 


(A-9b) 


(A- 10) 


where  is  in  atmospheres.  The  Sutherland  relation  for  the  viscosity  of  air 
is 


JL 

**1 


(T^/Z  + 198.6 

rf  ) T + 198.  6 


(A- 11) 


where  T is  in  degrees  Rankine.  Boundary- layer  paramenters  can  then  be 
found  from 


(A- 12a) 
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(A-  12c) 


where  the  dimensions  of  p^x  are  in  atmosphere-feet.  The  quantities 
q^,  and  are  found  from  Eqs.  (A-5)  and  (A-6).  Also 


Re 1/5  CH  = 0.0572 


fc) 


,3%^11/5 


(WTT) 


(A-  12d) 


In  a shock  tube  of  diameter  d,  operating  at  its  maximum  test  time,  the 
separation  between  the  contact  surface  and  the  shock  is  denoted  by  and  is 
found  from  (Ref.  2) 


1 *m  0.  1768  W (W  - 1) 


<P,d) 


W + 1.  25  W - 0.  80 


(P1X)1/5  fir) 


(A- 13) 


An  equivalent  flat  plate  Reynolds  number,  based  on  the  distance  that  a free- 
stream  particle  has  moved  relative  to  the  wall,  is**  "Re  = U£  (W  - 1)  x/ve- 
The  value  of  Re  at  the  contact  surface  is 


(P  4<*)' 


Re  = 6.  93  x 106  M — (W  - l)2  ^ L. 

m s u d . 

^e  (p.d) 


(A- 14) 


Turbulent  theory  is  applicable,  provided  is  of  order  10  or  greater. 


Re  is  not  to  be  confused  with  Re. 


APPENDIX  B 


APPROXIMATE  ANALYTIC  FORMULATION 


Approximate  analytic  expressions  for  boundary-layer  properties  can  be 
obtained  by  use  of  the  approximate  expressions  for  Eq.  (B-l)  presented  in 


For  ideal  gases,  the  latter  are 

0/6  . 3W  + 7 

1 - W 80 

Y = 7/5 

(B-la) 

. W + 2 
“ 24 

•< 

II 

un 

CO 

(B-  lb) 

^ n <C-7 

TT-w  - °*  157 

Y = 7/5,  W > 2 

(B-lc) 

= 0.  176 

Y = 5/3,  W > 2 

(B-  Id) 

Equations  (B-la)  through  (B- Id)  are  correct  to  within  8%  for  strong  shocks 
in  air  with  W < 15. 


Consider  the  case  of  a strong  shock  moving  into  air.  Neglecting  terms 
of  order  1/M^,  1/W,  and  h /h  compared  with  1,  we  find  h /h.  = 0.2  M^, 

S 6 6 1 9 

h /h  =1.9,  and  h /h  = 0.700.  Also,  for  T .=  522°  R and  198.  6/T  « 1, 

r e m e l . m 

the  Sutherland  relation  becomes  Pm/Vf  = 1*  380  (T^/T^)  ' . Then,  if  it  is 

assumed  that  T /T.  = h /h.  = 0.  14  M , the  viscosity  law  becomes  u /u ..  = 
m 1 m 1 s 7 m rl 

0.  524  M . Substitution  of  the  relations 


= 0.524  Ms 

h /h  = 0.  700 
m e 


(B-2) 
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and  Eq.  (B-la)  into  Eq.  (A- 12a)  yields 


(Pjx)1/5  i = 0.0370 


(W  - IK 


[W  + (7/3)]4 


1/5 


(B-3) 


Equation  (B-3)  is  in  remarkably  good  agreement  with  all  the  numerical 
results  in  Table  1,  including  W near  1,  since  6 /x  is  relatively  insensitive  to 
the  approximations  used  to  obtain  Eq.  (B-2).  Equation  (B-3)  agrees  with 
the  real  gas  results  in  Table  lb  to  within  5%  and  with  the  ideal  gas  results 
in  Table  la  (including  the  range  1.02  < W < 6)  to  with  10%.  Hence,  all 
boundary-layer  properties  of  interest  can  be  deduced  from  Eqs.  (B-l)  and 
(B-3)  with  an  accuracy  of  about  10%.  In  particular,  the  shear,  heat  transfer, 
and  mass  loss  to  the  boundary  layer  may  be  expressed  [from  Eqs.  (A-7) 
(B-l),  and  (B-3)] 


(Plx)1/5  Cf  = 2 Pr2/3  (Plx)1/5  C 


= 2. 22  x 10 


’r 


-3 
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(P,x)1//5  (-vg) 
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= 4.  65  X 10 


-3 


W + (7/3) 

L (W  - 1)ZJ 


(W  - 1)~ 


11/5 


(B-4a) 


[W  + (7/3)f 
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(B-4b) 


These  equations  agree  with  the  numerical  results  in  Table  1 to  within  about 

10%. 
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appendix  c 


WALL-SURFACE  TEMPERATURE 


In  Appendix  A,  it  was  assumed  that  the  wall- surface  temperature,  T^, 
remains  constant  at  its  original  value,  T^.  The  validity  of  this  assumption 
is  examined  herein  for  the  case  of  a turbulent  boundary  layer  behind  a 
moving  shock. 

The  transformation  t = x/u^  relates  the  variation  with  time  at  a fixed 
point  on  the  wall  (laboratory  coordinates)  with  the  variation  with  distance 
behind  the  shock,  x,  in  shock-fixed  coordinates.  Hence,  in  laboratory  co- 
ordinates, the  heat  transfer  to  the  wall,  as  a function  of  time  after  shock 
passage,  has  the  form  -qw~  t"1/5.  The  departure  of  wall-surface  tem- 
perature, T , from  its  initial  value,  T^,  can  be  found  from  (Ref.  6) 


T 

w 


(-yH175  ft  dr 
v/"V(p£k)b  JQ  (t— t)1^ 


(C-l) 


where  (pCk)b  is  the  product  of  density,  specific  heat,  and  thermal  conductivity 
of  the  wall  ("body")  material.  Dependence  of  the  latter  properties  on  tem- 
perature is  neglected.  Introduction  of  the  beta  function  8(0.5,  0.8)  = 1.297  /tT 
permits  Eq.  (C-l)  to  be  written  in  the  form 


1.297  (-qw)t1/2 

W"  Tl  = V^CkTb 


(C-2) 


Expressing  (-qw)  in  terms  of  Cpj  then  yields 
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(C-3) 


w 


Pl^x\Tl 


- 1 


1.297 

v/pCk)b 


*/m" 


(W 


Thermal  properties  of  typical  wall  materials  are  given  in  Table  2 in 
cgs-cal-°K  units.  In  similar  units 


= 0.05383 


TJT 


cal 

TIT 


cm 


sec 


atm9  K 


(C-4) 


for  air  at  a temperature  of  290°  K (522°  R).  Substitution  of  these  quantities 

into  Eq.  (C-3)  provides  the  variation  of  wall  temperature,  T , as  a function 

- 1 /5  . w 

of  distance  x in  cm.  Since  C„~  (p.x)  , it  follows  that  (T  - T.)~ 

0.8  0.3  HI  w 1 

Pj  x • 

In  the  limit  of  a strong  shock  in  air,  Eq.  (C-3)  can  be  expressed  in  the 

form 


Tnmrrfxf-  «)  = 


,)  = E M 5/2  W4/5  [ 

1 + o (4-) 

/ » l 

\W7 

(C-5) 


where  p,  and  x are  in  atmospheres  and  feet,  respectively,  and  E = 2.024  X 
-4  * 1/2 

10”  /(pCk),  ' . Equation  (C-5)  employs  Eq.  (B-4a)  and  assumes 

b 2 

h /h.  = 0.  38  M , h /h  « 1.  Values  of  E are  tabulated  in  Table  2 for 
r 1 s w r 4 

various  wall  materials.  The  quantity  E is  of  order  5 X 10  for  metal  walls 

_ 3 

and  of  order  5 X 10  for  nonmetal  walls.  The  boundary- layer  solution  in 

Appendix  A assumed  Tw  = Tj.  The  increase  in  wall-surface  temperature 

with  x has  a small  effect  on  this  boundary- layer  solution  if  Tw  - Tj  is  small 

relative  to  the  recovery  (adiabatic)  temperature,  T . Assuming  T /T.  = 

2 r r i 

0.  38  M , we  can  express  Eq.  (C-5)  by 
8 
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In  general,  (Tw  - T^J/T.  is  small;  thus,  the  variation  of  wall-surface 
temperature  need  not  be  considered  in  the  boundary- layer  equations.  How- 
ever, (Tw  - Tj)/Tj  may  become  sufficiently  large  that  melting,  or  ablation, 
occurs  for  very  strong  shocks.  In  the  latter  case,  Eqs.  (C-5)  and  (C-6)  are 
inapplicable.  For  example,  if  M = 10,  W = 8.  3,  the  right-hand  sides  of 
Eqs.  (C-5)  and  (C-6)  become  1.0  and  0.03,  respectively,  for  metallic  walls. 
The  present  solution  is  then  consistent  and  applicable  for  p^  of  order  1 atm 
and  x of  order  1 ft.  However,  for  nonmetallic  walls,  under  the  same  flow 
conditions,  the  right-hand  sides  of  Eqs.  (C-5)  and  C-6)  are  approximately 
10  and  0.  3,  respectively,  and  some  ablation  may  occur  at  distances  on  the 
order  of  1 ft  or  greater  behind  the  shock.  Ablation  will  tend  to  limit  the 
wall-surface  temperature  rise  and,  thus,  make  the  present  assumption 

h /h  « 1 valid,  but  the  effect  of  wall  blowing  must  be  taken  into  account, 
w r 

This  is  done  in  Appendix  D. 


WALL  BLOWING  EFFECTS 


The  effect  of  blowing  on  wall  shear  and  heat  transfer  is  not  well 
established  for  turbulent  boundary  layers.  No  experimental  data  are  avail- 
able for  the  case  of  a turbulent  boundary  layer  behind  a moving  shock.  Never 
theless,  an  estimate  of  this  effect,  based  on  results  for  semi-infinite  flat 
plates,  is  presented  herein.  The  blown  gas  is  assumed  to  be  air -like. 

8 9 

Based  on  the  flat-plate  analogy,  ' the  ratio  of  blowing  to  nonblowing 
(subscript  zero)  shear  and  heat  transfer  may  be  expressed  by 


W.  H.  Dorrance,  Viscous  Hypersonic  Flow,  McGraw-Hill  (1962) 
pp.  59,  206-220. 


W.  M.  Kays  and  R.  J.  Moffat,  "The  Behavior  of  Transpired  Turbulent 
Boundary  Lavers,  " in  Studies  in  Convection,  ed.  B.  E.  Launder, 
Academic  Press,  New  York  (1975),  p.25l. 


Note  that 


B s In  (1  + B') 

If  we  assume  that  all  the  thermal  conduction  from  the  air  to  the  wall 
results  in  ablation  of  wall  material  (i.  e.  , thermal  conduction  within  the  wall 
is  neglected),  it  follows  that 


= (pv)w  AH 


(D-3) 


where  AH  is  the  heat  of  ablation  of  the  wall  material.  (More  generally,  AH 
is  the  difference  in  enthalpy  between  the  unheated  wall  material  and  the  wall 
material  at  the  air/wall  interface. ) Equation  ( D-2b ) can  then  be  expressed 
as 


B' 


(D-4) 


and  the  effect  of  blowing  on  wall  shear  and  heat  transfer  can  be  readily 

deduced  from  Eq.  (D-lb).  The  value  of  B*  can  be  readily  estimated.  For 

air  at  an  initial  temperature  of  522°  R,  h.  = 69.4  cal/gr.  Thus,  for  strong 

2 2 1 

shocks,  hf  = 0.  38  Mg  h^  = 26.4  Mg  cal/gr.  For  a vaporizing  (T  = 3100®  K) 
steel  wall,  AH  = 2100  cal/gr  and  B'  = 1.25,  5.00,  and  125  for  Mg  = 10,  20, 
and  100,  respectively. 

The  variation  of  boundary-layer  momentum  thickness  with  distance  is 
of  interest.  The  momentum  integral  equation 


2 d9  i / v / \ 

p u 3 — = T + (pv)  (u  - U ) 

e e Ax  w r w e w 


(D- 5a) 


can  be  expressed  in  the  form 


- i de  _ f , vpv,w 
(W  - l)z  ^ ^ Peue  (w  - 


(D-  5b) 


Taking  the  Reynolds  analogy  to  be  = 2C^,  substituting  from  Eq.  (D-  lb),  and 
integrating,  we  obtain,  at  a fixed  x 


0 _ (1  + B')  In  (1  + B') 


(D-6) 


from  which  momentum  thickness  can  be  obtained.  The  vertical  velocity  at 
the  edge  of  the  boundary  layer  is  also  of  interest.  The  latter  can  be  expressed 
by 


% = dS* 
u dx  p u 


(D-7) 


5JC 

Assuming  is  independent  of  x and  recalling  d6Q/dx  = (4/5)  (6q/x), 

we  can  express  Eq.  (D-7)  as 


In  (1  4-  B / ) (plx)  *CH*0 

1 6*  (bQ  /b0\  (P1X)1/5  60 

'7*  \TT^/  - 


(D-8) 


The  ratio  6*/6q  is  needed.  If  it  is  assumed  that  the  form  factors  6*/6  and 
9 /«  are  unchanged  by  blowing  (which  is  highly  unlikely),  then,  at  a fixed  x 


6 _ 0 _ (1  + B')  In  (1  + B') 

r 


(D-9) 


The  validity  of  Eq.  (D-9)  requires  further  study.  In  the  absence  of  better 
estimates,  Eq.  (D-9)  can  be  used  to  evaluate  (,*/(,*  in  Eq.  (D-8).  Substitution 
of  Eqs.  (B-lc),  (B-3),  and  (B-4a)  into  Eq.  (D-8)  then  yields,  for  air  (Pr  = 0.9) 


e (1  + B')  In  (1  + B/) 


0.256B'  /w  + (7/3)\ 
’ 1 + B'  ^ W - 1 J 


= (l.-jL0.744B-nn(l+B-)  for  w >:>  t 


(D-lOa) 

(D-lOb) 


When  W»  1,  wall  blowing  increases  v /(v  )-• 

C 0 0 

The  rate  of  growth  of  the  boundary- layer  thickness,  6,  for  blowing  and 
nonblowing  cases  can  also  be  compared.  Let  x and  xQ  denote  the  correspond- 
ing axial  stations  for  blowing  and  nonblowing  cases,  respectively,  where  the 
boundary- layer  thickness  has  the  same  value.  Since  6 x4/5,  Eq.  (D-9) 

indicates  that  for  fixed  6 


B' 


7TT  ByTTn  ( 1 + B * ) 


5/4 


(D-ll) 


Equation  (D-ll)  is  useful  for  estimating  the  effect  of  blowing  on  boundary- 
layer  closure  in  tubes. 
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SYMBOLS 


B,B' 

crcH 


Re,  Re 


speed  of  sound 

wall  blowing  parameters,  Eqs.  (D-2a)  and  (D-2b) 

local  shear  and  heat  transfer  coefficients, 

Eqs.  (A-7a)  and  (A- 7b) 

static  enthalpy  per  unit  mass 

recovery  (adiabatic  wall)  value  of  h in  laboratory 
coordinates,  Eq.  (A-2c) 

shock  Mach  number,  u /a, 

w 1 

Prandtl  number 
pressure  in  atmospheres 
heat  transfer  to  wall,  (k  8T/dy)w 
Reynolds  numbers 

1 /3 

recovery  factor,  (Pr) 

temperature,  degrees  Rankine 

velocity  in  x direction,  shock-fixed  coordinates 

shock  speed  (wall  velocity  in  shock-fixed  coordinates) 

free- stream  velocity  behind  shock  in  shock- fixed  coordinates 

vertical  velocity  at  edge  of  boundary  layer 

uw/ue'  Pe/pi 
streamwise  distance 

ratio  specific  heats 

boundary- layer  thickness 


displacement  thickness,  shock-fixed  coordinates 


i 


0 momentum  thickness,  shock-fixed  coordinates 

P viscosity 

u kinematic  viscosity 

p density 

Tw  wall  shear,  (ji8u/ay)w 


Subscripts 

e 

m 

r 

w 

0 

1 


free  stream  downstream  of  shock 
mean  value 
recovery  value 

wall  value  downstream  of  shock 

zero  wall  blowing  value 

free- stream  value  upstream  of  shock 


THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structur.il  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
o:  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnctospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  spac  e systems. 
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